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The multipole response of ^Hbat drops doped with very attractive impurities, such as a Xe atom 
or an SFg molecule, has been investigated in the framework of the Finite Range Density Functional 
Theory and the Random Phase Approximation. We show that volume (L = 0) and surface [L = 1, 
2) modes become more fragmented, as compared with the results obtained for pure ^Hejv drops. In 
addition, the dipole mean energy goes smoothly to zero when A'^ increases, indicating that for large 
N values these impurities are delocalized in the bulk of the drop. 



The study of helium drops has been the object of extensive experimental and theoretical investigationsETa. One of 
the goals of these studies is to understand how various bulk physical properties of the quantum liquid are modified 
in restricted geometries. Special, mention deserves the emerging field of infrared spectroscopy of molecules inside 
or attached to helium dropletscTH which has motivated a W^&t theoretical activity to determine how the molecular 
moments of inertia are affected by the helium environmentQ'El. 

The main experimental effort has focused on the study of puro-|ajid-|doped ^Hcat drops, for which a microscopic 
description of the ground sta1je||fe) using Monte Carlo techniqueg3'Ert3, and of the elementary excitations using an 
optimized variational methodlijO are available. Density functional calculations of the gs and excitation spectrum 
using finite-range (FRDF) or zero-range density functionals have been carried out, see Refs. and references 

therein. Recently, the physical apaearance of quantized vortices pinned by dopant molecules in "Tie droplets has been 
studied within the FRDF theoryEj. 

In the case of '^Hc drops, experimental data are becoming availableEi, Small '^He drops are difficult to produce since 
a minimum number of atoms is needed to produce a selfbound dropcJ, and are as difficult to observe as *He drops 
because they are neutral. Nevertheless, ^He systems constitute the only Fermi systems capable of being observed 
in bulk liquid and droplets, and for this reason they have attracted some ptheoretical interest. Yet, microscopic 
calculations of '^He droplets are scarce, and only concern the structureLjl£3. A mass-formula for ^He drops 
based on an extended Thomas-Fermi method has been proposed — , and the binding energy of open-shell '^He drops 
has been determined by a semiempirical shell-correction methpdE^. The gs of small polarized Li-^HeA? clusters has 
been determined using the Path Integral Monte Carlo MethodEJ. Gropjd state properties of '^He drops doped with 
some inert atoms and molecular impurities have beeU|Xecently studied^ within the FRDF theory, as well as the gs 
structure of pure or doped mixed '^He-'*He dropletsHt^. Studies of mixed droplets are relevant in connection with 
the experimental results presented in Ref. ^. Indeed, it is crucial to know the composition of the first solvation shells 
around the impurity to determine if the dopant molecule is in a superfluid environments, or to determine whether the 
molecule may couple to bosonic or fermionic-type liquid excitations which in turn determines the dissipative picture 



Previously quoted references indicate that there has been an enormous impetus in the development and application 
of microscopic techniques to the description of liquid helium drops. However, current experiments sometimes have 
to deal with situations that cannot be addresaad-by fully microscopic methods. We can mention, for example, the 
description of very large '^He and ^He dropgE3c3, or the structure of large mixejd, drops already discussed. As a 
matter of fact, in spite of the recent progress made in the variational descriptionCJ of small (up to 40 atoms) ■^He 
droplets, a simultaneous description of ground state and elementary .excitations of pure '^He droplets has beeijunnly 
obtained within the density functional theory, using either zero-rangd23i2j or finite-range density functional£jc3. In 
these situations, density functional theory results provide a useful guide to the appearance of interesting physical 
phenomena obtained at the price of introducing some phenomenology. 

The aim of this work is to analyze the distortions caused by the presence of an impurity like a Xe atom or an SFg 
molecule in the excitation spectrum of a '^He droplet. The solvation energies of these impurities have been found to 
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be negativecll, and this makes plausible the scenario underlying in our calculation. We have analyzed this effect in the 
framework of the FRDF theory and the Random Phase Approximation (RPA). This paper is organized as follows: 
In Section II we briefly introduce the finite-range density functional we use, and the particle-hole (ph) interaction 
employed in the RPA calculations. In Section III we present results for the volume L — and the low-multipolarity 
surface excitations. A preliminary account of the dipole response has been previously reportedlj. In Section IV we 
draw the conclusions, and in an Appendix we present an example of how the angular decomposition of the ph matrix 
elements has been carried out. 



II. THE FINITE RANGE DENSITY FUNCTIONAL AND PARTICLE-HOLE INTERACTION 



In the framework of the density functional theory, the ground state of ■^He^r doped drops is found by minimizing 
the energy E written as 



wheri 
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E[p,t] = J dr£{p,T) , 
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The particle p{r), current j{r) and kinetic energy r(r') densities are written in terms of the single particle (sp) 
wave functions (pkir) obtained solving the Kohn-Sham (KS) equations deduced from Eq. (||). For systems having 
an effective mass m*{r), the inclusion of a term j'^/p in Eq. (H) guarantees that the density functional is Galilean 
invarianto. This term has no influence on the ground state of time-reversal invariant, spin-saturated droplets, and 
for this reason it is usually omitted. However, its contribution to the ph interaction for systems like the present one, 
in which the impurity is treated as an external field breaking the translational invariancc of the system, cannot be 
neglected. 

In the above expression VLj{\r — f'\) is the Lennard- Jones interatomic potential screened at short distances 
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and the averaged density p{f) is defined as 
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The effective mass m* is parametrized as m* = rn(l — p/pc)~^- The set of coefficients entering the definition of £{p, r) 

in the case of Xe, and 
In both cases we have assumed that the impurity 



can be found in Table I of Ref. |3l|. Vimp is the helium-impurity potential taken from Ref. 
from Ref. ^ in the case of SFg, in its spherically averaged version, 
is an object of infinite mass located at the coordinate origin. 

The distortion of the ground state structure of "^He drops due to the presence of impurities has been described in 
detail in Ref. ^ To analyze the multipole excitations induced by an external field that couples to the particle density 
of the drop, we have used the time-dependent version of the density functional theory. For sufficiently weak external 
fields the response can be treated linearly within the RPA. In this approximation the elementary excitations of the 
system are described in terms of correlated ph transitions. The amplitude of a particular excited state in the basis 
of a discrete space of ph transitions is obtained by diagonalizing the Hamiltonian H = Hq +,V„h. which is the sum of 
the KS Hamiltonian Hq plus the ph interaction Vph ■ This is done solving the RPA equatior 
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where the matrices A and B are written in terms of matrix elements of the interaction between the ph pairs that can 
be coupled to have the desired angular momentum. 

Writing the particle, current and kinetic energy densities in terms of the sp basis (j)k{'P) and the occupation numbers 
Pki 

P{r) ^^(l)l{r) Pki (1)1 (r ) (7) 

kl 

= ^ - ^U^')Pkl Mr) \r=r' (8) 

kl 

r{r ) = ^l^^)Pki Ur')\r=r' , (9) 

kl 

the ph interaction is obtainedSEl from the second variation of the energy functional with respect to the occupation 
numbers: 

V.jki = {ij\Vph{ri,r2)\kl) = ^ f . (10) 

OPtkOPjl 

If m* — m in the density functional (a situation which we indicate with the notation E — E[p\), the second variation 
of the energy with respect to the occupation numbers taken at the ground state straightforwardly provides the ph 
interaction Vph- This variation can be obtained as 
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and comparing with Eq. (HO) it results 
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The presence of a position-dependent effective mass in the functional introduces a velocity dependence in the ph 
interaction. In this case Eq. (|l^) becomes 
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The arrow on the gradient operators indicates whether they act on the left or on the right. The terms arising from 
current derivatives are essential to fulfill the Thomas-Reiche-Kuhn (or energy-weighted) sum rule. The contribution 
of current terms to the ph interaction is 
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where a sum over the three components a is assumed, and the gradients only act on the sp wave functions. This 
expression coincides with the back-flow contribution to the ph interaction for ^He dropsu. 
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accounting for Eqs. (p^, ( p^ ) and that 



Particularizing to the density functional Eq. 
( |l3| ) gives for the ph interaction 
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Equation ( [l6| ) shows that in addition to the Lennard- Jones potential Vlj, Eq. (^, the ph interaction has finite-range 
terms, velocity dependent components, and other terms which combine both finite-range and velocity dependence 
through the presence of gradient operators. 

The next task is to calculate the matrix elements in the ph basis. This is greatly simplified in the case of droplets 
with a magic number of ^He atoms, the only droplets studied here. In this case, the mean field is spherically symmetric 
and the angular part of the sp wave functions is a spherical harmonic. Performing a multipole expansion of the ph 
interaction, the sum over third components can be done and only radial integrals remain to be numerically computed 
(see the Appendix for details). This allows one to compute the RPA matrices A and B. After diagonalizing Eq. (H), 
the strength function from the gs |0) to the set of excited states {|n)} (with excitation energies {wn}) is obtained as 



MQLm 



(19) 



where is the excitation operator for which we have made the natural choices Qq — rf for L — Q (volume mode), 
and Ql — Ti^Loi^i) for L = 1, 2 (surface modes). 

The transition matrix element of the strength function is obtained in terms of the solutions of RPA equation and 
the sp radial wave functions Ui defined in the Appendix. The explicit expression for surface modes is given by 
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where (i'm||yL|Kj) is the angular reduced matrix element of the excitation operatoiEl. The corresponding expression 
for the monopole mode is 

(OlQoln) = ^(X^^) - y!^;>) {ura\r^\u,) ^t^u ■ (21) 

The transition (also called induced) densities for the operator Ql — X^i ''f ^io(f^i) that causes surface excitations 
are obtained as 

p„o(0 ^ --=^E(-^^^ - ^i."^) ^U\yL\%) ""('^^^f , (22) 

rai 

and (OlQ^jn) = J dr r^~^^ pno{r). The corresponding induced densities for the monopole mode are given by 

Pnoir) = 2^(a;„^ - y„\/) bi^i^ , (23) 

mi 

and {0\Qo\n) = J drr^pno{r). 

Obviously, the dimension of matrices A and B depends on how many particle-hole pairs mi are taken after dis- 
cretizing the continuum. We have included enough sp states so that the Thomas-Reiche-Kuhn sum rule is satisfied 
within 98%. We have also checked that for pure drops the dipole mode is at zero energy due to the translational 
invariance of the system. 

To finish this Section, we would like to recall that originally, density functionals for liquid "^He were obtained from 
a contact, velocity-dependent '^He-^He effective interactionEll that made it rather simple to evaluate the contribution 
of direct and exchange terms to the total energy and to the ph interaction. Later on, a finite-range component was 
added to the contact interaction to improve its properties at finite momentum. This is the origin of the screened 
Lennard-Jones potentialE3o, which takes care of two major characteristics of the interatomic potential the original 
effective He-He interaction lacked, namely the hard core repulsion at short distances, and the asymptotic behavior. 
Thus, exchange effects, which are known to be large in liquid ^He, are phenomenologically accounted for in the density 
functional through the effective parameters entering its definition. 



III. RESULTS 



A. Monopole mode 

Figure |l| shows a comparison between the monopole ('breathing mode') spectrum of pure and doped drops. It is 
seen that the presence of the impurity increases the fragmentation of the spectra in the high energy region. This 
effect is more important for small clusters and more attractive impurities. In both cases of pure and doped drops, the 
mean energy defined as u) = lj jS{ujj)/ S{ujj) lies above the atom emission threshold (Fermi energy changed of 
sign) and decreases as the number of ^He atoms of the drop increases (see Figs. || and ||). It is worthwhile to recall 
that for jpwjKe, ijJe jv droplets, except for rather small N values the monopole strength is in the discrete region of the 
spectrumll3't3'E3, and that the presence of a Xe or SFg impurity, also increases the fragmentation of the spectrum; for 
small drops the monopole strength lies in the continuum regionE3E3. 

We display in Fig. |^ the transition densities corresponding to the more intense monopole peaks of Xe-|-^He4o. 
This figure shows the well-known fact that the monopole is a volume mode: the induced densities have a node and 
penetrate inside the drop. The bulk oscillations are connected with the oscillations in the drop density p(r), also 
shown in Fig. |[ which are due to the distribution of '^He atoms in solvation shells around the dopant on the one 
hand, and to the repulsive core of the effective interaction, on the other hand. 



B. Dipole mode 



The L = 1 spectrum shows again that fragmentation increases for the more attractive impurities (see Fig. ^). In 
this case the mean dipole energy always lies below the continuum threshold and decreases with the number of "^He 
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atoms of the drop. In spite that small doped drops are stable to dipole fluctuations since a large energy is needed 
to induce the oscillation of the impurity against the '^He atoms, Fig. |3| shows that when the drop size increases the 
dipole mean energy rapidly decreases and the mode eventually becomes unstable. This is considered a clear signatuxa 
that the impurity is delocalized in the bulk of the dropOScj. The dipole mode has also been found to be unstableEHl 
for large "'Hejv drops doped with inert atoms and SFg, for which the mean dipole energy lies in the discrete part of 
the spectrum. 

C. Quadrupole mode 

Fig. 1^ shows that, as compared to the pure case, the presence of a rather attractive impurity pushes this mode 
downwards in energy. When this causes the quadrupole mode to move from the continuum to the discrete part of the 
energy spectrum, the fragmentation decreases and the peak becomes more collective. This is the case for N — 40, for 
example. For larger drops, the quadrupole mode is below the atom emission threshold, see Fig. ^ (this also happens 
in pure drops), and the effect is not so clearly seen. 

Examples of induced densities for dipole and quadrupole modes are shown in Fig. 0. They are localized at the drop 
surface, as it corresponds to the surface character of these modes. 

We would like to close this Section indicating that in the case of pure "^He droplets, a cxumiarison with results for 
L = and 2 modes obtained using density functional built using fairly different strategiesH'Ej yields an overall good 
agreement. 

IV. CONCLUSIONS 

We have investigated the multipole collective excitations of '^Hejv drops doped with Xe atoms and SFg molecules 
in the framework of the FRDF theory plus the RPA. A comparison with the results for pure drops shows that the 
presence of these strongly attractive impurities increases the spectrum fragmentation. This effect appears in volume 
and surface modes as well, and it is more marked for small clusters and more attractive impurities. 

The m£sence of an attractive impurity decreases the mean energy of surface modes as in the case of doped ■^Hcat 
clustersE3. For large clusters the mean energy of surface modes lies below the atom emission threshold, whereas for 
the monopole volume mode it is always above the threshold. 

When the cluster size increases the dipole mean excitation goes to zero, indicating that the impurity is delocalized 
in the bulk of the drop for '^He clusters doped with Xe and SFg impurities. A similar effect was found in *He clusters. 
From the experience gathered in the case of *He clusters, we may conclude that whereas the precise value of the 
(rather fragmented) collective modes may be sensitive to the arbitrariness introduced in the choice of some of the 
FRDF ingredients, as for example the core of the screened Lennard- Jones potential, we consider robust the prediction 
of the impurity delocalization, as well as the evolution of the mean mode energies with the number of atoms. 
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APPENDIX A: 

For a spherically symmetric system the dimension of the ph space can be drastically reduced by analytically summing 
over the degenerate third components of the angular momentum. We illustrate this point taking as an example the 
Lennard-Jones contribution VLj{r) to the ph interaction assuming that the ph states are coupled to yield an orbital 
angular momentum {L, Ml), and a spin {S, Ms). Using the appropriate sp quantum numbers we represent the orbital 
(j)a in coordinate and spin spaces as 

(pa = (l)aeat^Ar)x:La^ ■ (Al) 
The matrix elements of the residual interaction (llq) between the ph states can then be obtained as 
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{mj\Vph\m) = ^(-1)^—'+^--^ - a,\SMs) (~c7„ - aj\SMs) {x.Jx.,) (x^Jx^J x 

cr's 

(£„^^/,„^ _ f^,\LML) {Ujt^n - t^j\LML) X 

/i's 

Due to the spherical symmetry of the mean field, the sp wave functions separate in radial and angular components: 

Mr) = ^^Y,^{f) , (A3) 
and one can perform a multipole expansion of the ph interaction 

Vph{ri2) ^Y.^L{rur2)Yl,,{h)YLM{f2) ■ (A4) 

LM 

The expression for the Lennard-Jones term of Eq. (|l|) can be finally written in terms of the reduced matrix 
elements of the spherical harmonica^ as 

{mj\VLj\m) = 25so^=^ {£m\\YLm) (^111111^,) , (A5) 

where /^^^^j is the radial integral 

nOO POD 

-^mjm = / dri / dr2U„ii^{ri)uje^{r2)VL{ri,r2)uni^{r2)uii,{ri) . (A6) 
Jo Jo 

These integrals are obtained numerically. To describe the states above the continuum threshold we have followed the 
usual prescription of enclosing the system in a sphere of large radius and require that the radial wave functions vanish 
at this distance. In this way we obtain a discrete spectrum of states that replaces the continuum. We have checked 
that the results are stable against reasonable changes of the radius of the sphere, which we have taken to be about 3 
times larger than the mean square radius of the cluster. 
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FIG. 1. Comparison between the monopole spectrum of pure and doped drops. Each excited state is represented by a 
vertical stick whose height gives its fractional contribution to the energy weighted sum rule. The arrows indicate the position 
of the atom emission threshold in each case. 
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FIG. 2. Same as Fig. |l| for the monopole spectrum of drops doped with Xe. 
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FIG. 3. Mean excitation energies oj (K) and chemical potential changed of sign ep as a function of A'^ for ^HeAr drops doped 
with Xe and SFg. The lines have been drawn to guide the eye. 
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Fig. 4 

FIG. 4. Transition densities (arbitrary scale) corresponding to the more intense monopole states of Xc+^He4o. The ground 
state density p(r) is also shown. The transition densities have been scaled to have a common value at r = 10 A. 
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FIG. 5. Same as Fig. nl for the dipole spectrum of doped drops. 



13 







Xe@3He24o 











0.4 



0.2 



0.0 





SFg@3He4o 




' i 


n-i-f [ 





6^ ' '^240 



0.2 



0.0 



2 4 

6;(K) 



1 2 

co{K) 



6 
Fig. 6 

FIG. 6. Same as Fig. ^ for the quadrupole spectrum. 



14 




Fig. 7 



FIG. 7. Transition densities (arbitrary scale) corresponding to the more intense L = 1 and 2 peaks and ground state density 
p(r) of the '^He4o drop doped with Xe and SFe. Scaling factors as in Fig. have been used. 
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